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Abstract. We present Spitzer 8 /im transit observations of the extrasolar planet system HD 149026. 
At this wavelength, transit light curves are weakly affected by stellar limb-darkening, allowing 
for a simpler and more accurate determination of planetary parameters. We measure a planet- 
star radius ratio of R p /R t =0.05158±0. 00077, and in combination with ground-based data 
and independent constraints on the stellar mass and radius, we derive an orbital inclination of 

■ _■ o 

£*H' i =85? 4 +0 o 9 and a planet radius of 0.755 ± 0.040 Rj. These measurements further support 
Q . models in which the planet is greatly enriched in heavy elements. 
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1. Introduction 

oo : 

t— ( , The hot, Saturn-mass exoplanct HD 149026b (M — 114M ffi ; discovered by Sato et al. 

2005) has garnered much attention due to its potential to directly test models of planet 
formation. The planet's small observed radius for its mass implies that roughly 2/3 of its 
mass is in the form of heavy elements- more than found in all of the Solar System planets 
combined. The existence of such a metal-laden planet orbiting a very metal-rich host star 
([Fe/H] = 0.36) is a smoking gun for the core-accretion theory (e.g., Pollack et al. 1996), 
though the immense quantity of heavy elements is still a challenge to canonical models of 
core accretion. HD 149026b has also received notice for its scorching day-side brightness 
temperature (2300 K), which is well in excess of its expected equilibrium temperature. 
The planet may be a standard bearer for the proposed 'pM' class of planets (Fortney et 
al. 2008), which are characterized by hot stratospheres and large day/night temperature 
contrasts caused by opacity at visible wavelengths due to gaseous TiO and VO high in 
their atmospheres. 

Unfortunately, the present fractional uncertainty in the key observable parameter, the 
planetary radius R p , is 7% (Winn et al. 2008). This study is inspired by the potential of 
infrared photometry with the Spitzer Space Telescope to reduce this uncertainty. Because 
of the near absence of stellar-limb darkening in the infrared, transit light curve modeling 
is greatly simplified and gives results largely independent of assumptions about limb- 
darkening coefficients. We present Spitzer 8 (im observations of the transit of HD 149026. 
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Together with previously published ground-based data, we derive precise constraints on 
the radius of HD 149026b and other system parameters^ 



2. Observations and Analysis 

We observed the UT August 14, 2007 transit of HD 149026 through the IRAC 8/im 
channel on the Spitzer space telescope (see Figure 1). We obtained 67,008 images in 
IRAC subarray mode at a 0.4 s cadence. The data exhibits very low levels of time- 
correlated, or "red" , noise, as demonstrated by Figure 2. Figure 2 also indicates that, at 
15 minute resolution (bins of roughly 2,000 images), the photometric precision dips below 
2 x 10~ 4 . In order to achieve the strongest possible constraints on the system parameters, 
we combined our data with previously published ground-based photometry (Sato et al. 
2005, Charbonneau et al. 2006, Winn et al. 2008). 

We employed Markov Chain Monte Carlo to derive best-fit parameters and uncertain- 
ties (see Table 1). We imposed an external constraint on the stellar mass, as derived by 
Sato et al. (2005). We also employed the empirical V, K relations of Kervella et al. (2004) 
for the angular diameters of dwarf stars, and combined with the Hipparcos parallax to 
derive an external constraint on the radius of the star HD 149026. 

IRAC 8/im observations exhibit a detector ramp (see Figure 1 and e.g., Knutson et 
al 2007). Our MCMC model includes a simultaneous analysis of a ramp model and the 

f For a detailed description of our analysis, please consult Nutzman et al. (2008). 
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Figure 1. Spitzer transit photometry for HD 149026, with 40 second resolution (bins of 100 
images). The top panel displays the raw light curve and exhibits the well-known IRAC 8 fim 
detector ramp (e.g., Knutson et al 2007). The middle panel displays the ramp-corrected light 
curve together with the best-fitting transit model. At bottom are the residuals from the best-fit. 
The root-mean-square residual is only 15 % greater than the expected photon noise. 
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Estimates of the HD 149026 
System Parameters 



Parameter Median 



15.9'" Perc. 



84. 1 th Perc. 



Rp/R* 


0.05147 


-0.00077 


+0.00076 


i [deg] 


85.3 


-0.8 


+0.9 


a/R* 


6.20 


-0.25 


+0.28 


R P [R.i] 


0.755 


-0.040 


+0.040 


R*\R®] 


1.497 


-0.069 


+0.069 



Table 1. Best-fit parameters and uncertainties of HD 149026 system properties derived 
through a Markov Chain Monte Carlo analysis. 

transit model, hence our error bars include the resulting correlations and uncertainties. 
The data, with and without the ramp correction and binned 100:1 are shown in Figure 
1, together with the best-fitting transit model. 



3. Discussion 

We have presented and analyzed Spitzer 8 /jm transit observations of the HD 149026 
system. By incorporating previously published data, and adopting constraints on the 
stellar mass and radius, we improve the determination of the planetary radius to R p = 
0.755 + 0.040 Rj. 

Our measurement reinforces previous findings of the intriguingly small radius of HD 
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Figure 2. The root-mean-square of binned residuals vs. bin size. The solid line is proportional 

to iV- 1 / 2 . 
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Figure 3. The radius of HD 149026b for three planet models of Fortney et al. (2006). The 
light grey band represents the observational radius constraint from this work. The dark grey 
box includes the age constraint of Sato et al. (2005). The "icy" model is for an H/He envelope 
and a 77Afe core of ice, while the "rocky" model is for an He/He envelope and a 65-M® core 
of olivine. A solar composition model (dot-dashed line), in which the planet is predominantly 
H/He, is excluded by the observational constraints. 

149026b. In comparison, models in which HD 149026b is composed purely of H/Hc require 
a radius greater than 1.1 Rj (see Figure 3). 
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